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1 Introduction 

Recently, the rapidity distributions of the charged particles and the ratios of 
the numbers of anti-hadrons to hadrons in nucleus-nucleus collisions have been 
investigated at SPS[1,2] and RHIC[3]. Also, the possibility of the formation of a 
quark-gluon plasma state has been studied very intensively. 

The nucleon-nucleus(7V-^) and nucleus-nucleus(^-^) collisions have been stud- 
ied from the multiple-scattering view[4,5,6]. The basic view of the multi-chain 
model(MCM) was discussed in [7]. Kinoshita etal. investigated this view in detail 
and succeeded phenomenologically to reproduce the old data on the multiplicity 
and the inclusive spectra of the leading nucleon and secondary particles etc. in 
cohisions of 7V-A[8] and A-A[9,10]. In MCM[11], it is assumed that the nucleon- 
nucleon(A'^-A'^) interaction is considered as a exchange of one multi peripheral chain 
from which secondary hadrons are emitted as shown in Fig.l. However, it was dif- 
ficult to perform analytical calculation for A- A collisions because of the complexity 
of the general formula in resolving the full combination of chain configuration. Also, 
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Fig. 1: MCM for A-B collision with A = 2 and B = 3. Wavy lines 
represents the inelastic N-N interaction, namely, N+N ^ A^^+A^+hadrons. 

the probabilistic approach [12] and the leading cluster cascade model with a recur- 
rence equation[13] were studied for N-A collisions. The essential point of MCM 
with successive collisions is that the number of N-N collision is equal to the mass 
number of nucleus A. Namely, nucleus looks as if nucleons in nucleus stand in a 
line. 

In previous paper [14], we investigated the extended MCM(EMCM) on the ba- 
sis of MCM with successive collision[15] and estimated the full combinations of 
chain configuration in terms of the vector-operator method in the moment space by 
consulting the Mellin transformation[13]. We derived the analytic forms for single- 
particle distributions in A-A collisions. This approach may apply to the N-A, but 
does not extend to A-A collisions. 

In this paper, we propose the new EMCM, namely, a unified treatment of N-A, 
N-A, A-A and A-A collisions. We introduce the concept of the operator matrix in 
the moment space to characterize these collisions. By means of the operator matrix, 
we calculate the single- particle distribution of the inclusive process A-\-B C-\-X. 
We assume the dynamics to scale with energy. A nucleus is treated as a set of 
mutually independent nucleons. The cascading of the produced hadrons is neglected 
owing to the long formation length in nucleus. Also, the longitudinal motion is 
treated. The transverse momentum distribution is assumed to be independent of 
incident energy and longitudinal momentum. 

In Section 2, the EMCM are investigated and the vector-operator formalism 
in the moment space is given. The operator matrix is introduced to treat the 
A-A collision and the calculation method is given. In Section 3, the structure of 
N-N collision is studied. The single-particle distributions of N-A collisions are 
investigated and their analytic forms are derived. In Section 4, A-A collisions are 
studied. In Section 5, conclusion and discussion are given. The definition and 
properties of operator matrix are given in Appendix. 

2 Extended multi-chain model 

2.1 Moment space and vector-operator notation 

We consider A-B collisions. The mass numbers of the projectile nucleus A and 
target nucleus B are A and B, respectively. We pay attention to the single-particle 



distribution of the inclusive process A + B ^ C + X{pq{x)) defined by 

^,,^^^^J_aiA + B_^C + X)^ (1) 

where C is the produced hadron and x denotes the longitudinal momentum frac- 
tion(Feynman variable). We express the projectile(target) fragmentation regions 
as X > ( X < ). 

Wc assume that a nucleus is treated as a set of mutually independent nucleons. 
By the Mellin transformation of Eq.(l), the moment p^^(a) is defined by 



(a) = 
Jo 



p^'^ia) = I dxx"-'p^^{x). (2) 

The inverse Mellin transformation is defined by (x) = J^^^^ dax^~°' {a). 
Here, we employ a vector-operator notation[15]. The operator p^^{a) in the mo- 
ment space is defined by 

<qp^^(a)|iV>=p^«(a), 

where < C\ denotes the produced hadrons in the final state and |A'' > denotes the 
nucleon state in the projectile nucleus A or the target nucleus B. 

2.2 A-B collisions and operator matrix ^^^(6, a) 

The A-B collision is characterized by the N-N collisions between the nucleon 
in nucleus A and the nucleon in nucleus B. The peripheral reaction N + N ^ 
A^+iV+hadrons is considered where N in the final state denotes the leading nucleon. 
In the MCM with successive collisions, the nucleon in nucleus A successively collides 
with the nucleons in nucleus B.Then, the coUison number of times is limited to B 
and vice versa. Thus, the total collision number in A-B collision is AB. 

We treat the A-B collision at the impact parameter b and in the moment space 
a. We introduce the operator Qij(b,a) to denote the probability in the impact 
parameter b and the longitudinal momentum fraction distribution in the moment 
space of the collision between the i-th nucleon in A and the j-th nucleon in B where 
i = 1, 2, . . . , j4 and j = 1,2, . . . , B. The operator matrix Qj^g{h, a) to characterize 
the total structure of A-B collisions, is defined by the rectangular or square array 
of the operators Qij (6, a) with A rows and B columns 



QAB{b,OL) 



lib, a) Qi2{b,a) ■■■ QiB{b,a) 



QAi{b,a) QA2{b,a) ■■■ QAB{b,a) 



(3) 



By the means of the calculus of operator matrix given in Appendix, we may derive 
the informations on the single-particle distributions of ^ + i? C + A by the 
sum decomposition rules defined in Appendix from QABi^i <^)- Namely, we use the 
sum of the products of the row elements(Spr) and the sum of the products of the 
column elements(Spc) for a) in regions of x > and x < 0, respectively. 

The calculus of operator matrix gives the convenience of calculating the probability 
of one event in A-B collision. For example, this operator matrix of the A-B collision 
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Fig. 2: EMCM for A-B collision with ^ = 3 and B = 4. The intersection 
points represent A^-A'^ collision(«). The solid lines denote nucleon(A). 



with A = 3 and -B = 4 is shown by the diagram in Fig. 2 in stead of Fig.l. The 
intersection point in the diagram denotes the N-N collision between the nucleon in 
A and the nucleon in B. It is assumed that the time of the i-th nucleon in A and 
the j-th evolve from QiB{b, a) to Qii{b, a) and from QAj{b, a) to Qij(b, a) in Fig.2, 
respectively. The ordering of the operator Qij (b, a) obeys this time-evolution. 



2.3 Single-particle distribution in A-B collisions 

We consider the inclusive reaction A + B ^ N + X where N denotes the leading 
nucleon. We introduce the operator Qo{b,a) defined by 



Qo{b,a) = r^{b)Gia) + X{b)Jia). 



(4) 



for N-N collisions in nucleus. Here, X{b) = a^^i JdsJ dtTA{s)TB{i)S^{b + s - t) 
where cr^^ is the inelastic cross section of N-N collisions. Also, Ta{s) and Tsii) are 
the normalized nuclear thickness functions of nucleus A and nucleus B, respectively. 
The X{b) denotes the inelastic interaction probability while the 77(6) expresses the 
passing-through probability(r?(6) + A(6) = 1 owing to the probability conservation). 
Operators G{a) and J{a) characterize the momentum fraction distributions as 
defined [G{a), J{a)] = and G{a) = 1. The operator J{a.) expresses the inelastic 
interaction of A-Acolhsions. By considering both operators G{a) and J{a), we 
may calculate full combinations of the chain configurations with ease. In fact, it 
corresponds to MCM with two kinds of chain. 

In this case, the elements of the operator matrix with A rows and B columns 
are given by Qij{b, a) = Qo{b, a) where i = 1,2,...,^ and j = 1,2,. . .,B, and 
hence 



QABib,a) 



Qo{b,a) Qo{b,a) ■■■ Qo{b,a) 



_Qa(b,a) Qa(b,a) ■■■ Qo{b,a) 



(5) 



We pay attention to the projectile fragmentation regions(a; > 0) and thus sum 
the products of the row elements of Eq.(5). ^From Eq.(56) in Appendix, we get 



SprQAB(^a)=A(r7(6) + A(5)) 



AB-B 



{Qo{b,a)y 



(6) 
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Furthermore, we use the relation [G{a), j{a)] = 0. Eq.(6) reduces to 

AB-B , s 



k=0 



X 
1=0 



AB 



= AY^ Pm{AB; b) J2 m AB, m, B){G{a)f-\j{a))\ (7) 

m=0 i=0 

where Pm{AB;h) is the Glauber probability and H{l;AB,m,B) is the hypergeo- 
metric distribution. They are defined by 

/ AB — m\fm 

Pm{AB-h)=(^^\r^{b)^''-^\{hr, Hil;AB,m,B)=^ 



B 

^From the Glauber probability, we obtain the inelastic cross section of A-B collision 
given by 



AB 
inel 



„ AB 

SY, Pm{AB; b)= db[l - (1 - A(6))^«] 

m=l 



and the averaged collision number n = ABa^J^^ I'^tJ^i ■ 

The matrix elements of the leading nucleon N in the final state, is defined by 



1 



< N\J{a)\N >=^< N\p'''\a)\N >= F{a), (8) 

where < and > denote the leading nucleon state and the nucleon state in the 
initial nucleus A or B. By the inverse Mellin transformation, the fraction functions 
F{x) is given by 

1 />r+ioo 

F{x) = — / daa;i-"F(a), (9) 

where F{x) is normalized to unity, namely, ^F{x) = 1. 

The distribution Qj;jJj^(^,a,A) in the projectile fragmentation regions(A side) 
is defined by 

Qj,%{b,a,A) =< N\SprQMb,a)\N > . (10) 
When we put G{a) = 1, from Eqs.(7), (8) and (10), we get 

AB B 

Qj,%{h, a,A) = AY, Pm{AB; b) ^ H{1; AB, m, B)F{ay. (11) 
m=0 ;=o 
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By means of the inverse Melhn transformation of Eq.(ll) except for the passing- 
through terms, the single-particle distribution of A + B ^ N + X for the projectile 
fragmentation regions (x > 0) is given by 

/AB B 
dbA Pm{AB; b) H{1; AB, m, S)i^(') {x) (12) 

TO=1 (=1 

with 

F«(a;)= /'^i^(-)F('-i)(y) with FW(y) = (5(1- y). (13) 

Jx y y 

Eq.(12) agrees with the result discussed in [14]. From Eq.(12), wc get the inelastic 
cross section given by a^f; = ^ I ^Pn^{x) = I db[l - (1 - X{b))^^]. Similarly, 
we may investigate the N-A collisions and get the same result as the one in [14]. 



3 Single-particle distributions in N-A collisions 

3.1 Structure of N-N collision in nucleus 

We consider the N-N collision in nucleus. In this case, we must introduce 
the annihilation interaction between N and N. Then, the leading mcson(ML) is 
produced by the annihilation interaction as N + N ^ Mi + M^+hadrons. Thus, we 
newly consider Ml-N{N) collisions and M^-Ml collisions in nucleus. We consider 
the peripheral reactions Ml + N{N) Ml + iV(iV) -I- hadrons and Ml + Ml ^ 
Mi-l-Mi+hadrons, but neglect the non-peripheral reaction Ml+N{N) -> N{N) + 
Mi -I- hadrons because of its small probability. We introduce the following operator 
set: 

(i) N-N coUision 

(a) Qn{b,ct) = r]n{b)Gn{ct) + A„(6)J„(q;) for non-annihilation collision, 

(b) Xa{b)A{a) for annihilation interaction, 

where r;(6)„ + A(6)„ + A„(6) = 1. Also, J dbX^ib) = atn and / fi/7A„(?7) = af/ 
where tr^^ is the non-annihilation cross section and is the annihilation 
cross section in N-N collision. 

(ii) Ml-N{N) coUision 

(c) QL{b,a) = rjL{b)GL{a) + XL{b)JL{a) for inelastic collision, 
where r;(6)i + A(6)l = 1 and / dbXL{b) = af^^f^^^ 

(ii) Ml-Ml collision 

(d) Mo{b,a) = ri,n{b)Gm{a) + Xm{b)Jm{ot) for inelastic collision, 
where r;(6)„ + A(&)„ = 1 and / d5A„ (6) = a^'^i^\ 

The operators G„(a), GL(a) and Gm{ct) denote the passing-through. The 
annihilation operator A{a) is not commutablc with Qn{b,a) and QL(b,a). The 
operator ordering of A{a) is in front of (5l(6, a) and at the back of Qn{b,a). 
Notice that [G„(q;), J„(a)] = 0, [GL{a), ^/.(a)] = and [Gm{o:), Jm{a)] = 0. Also, 
we define 

< N\Jn{a)\N >=< N\Jr,{a)\N >= F{a), 

< ML\A{a)\N >=< ML\A{a)\N >= ^(0;), 
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Fig. 3: Diagram for the annihilation term in N-A coUision. The intersection 
points represent M^-N conision(n), N-N conision(») and annihilation in- 
teraction of N-N collision(o). The solid and dashed lines denote nucleon(iV) 
and leading meson(Mi), respectively. 

< N\JL{a)\N >=< N\JL{a)\N >= N{a), 

< ML\JL{a)\ML >= L{a), < ML\Jm{a)\ML >= M„(a) 

and the others of matrix element are zero. These functions in x space are normalized 
to unity as F{x). 

3.2 N-A collisions 

We investigate N-A collisions at the impact parameter b. We treat the inclusive 
processes oiN + A—i'N + X and Ml + X for the projectile fragmentation regions 
of N{x > 0) and the inclusive processes ofN + A^N + X and Ml + X for the 
target fragmentation regions of A{x < 0). The N-N and Ml-N{N) collisions in 
the operator set contribute to N-A collisions. The operator matrix Qjva(^j with 
one row and A columns to characterize N-A collisions at the impact parameter b 
and in the moment space, is given by 

QNAib, a) = Q%A{b, a) + a), (14) 

where QffAi^, a) and Ojva(^) denote the non-annihilation and annihilation terms, 
respectively. They are given by 

q!va(^ a) = [Qn{b, a), • • • , Q„(6)], (15) 
^ ' 



A-l 



QNAib, a) = A„(6)[Qi(6, a), ■ • • , Qi(6, a), A{a), Qn{h, a), • • ■ , Qn{h, a)]. (16) 

I A-l-l 

The diagram for the annihilation term in N-A collision is shown in Fig. 3. 

(1) N+A N+X and Ml+X for the projectile fragmentation regions(a; > 0) 



7 



We apply the sum decomposition rule of Eq.(52) in Appendix to Qna{^,oi)- 
^From Eq.(15), we have 

A 

SprQlv^(6,a) = {Qn{b,a) f = ^ &)(J„(a))'", (17) 

m=0 

where P„(A; b) ^ ( ^] r/„(6)'^-"A„(&)™ and then G{a) = 1. 



m 

From Eq.(16), we have 



A-l 



SprQ^^(6,a) = ^ A„(6)(Qi(6, a))'i(a)(Q„(6, a))^"'"! 



A-l I 



Ub)P^{l-APt{A-l-lMJL{a)yA{a){Ua))\ (18) 



(=0 s=0 t=0 



where P/'(/;6) = (g) ??L(b)'"'Az,(6)^ and then ^1,(0!) = G„(a) = 1. 
Taking the < -/V| . . . |^ > matrix element of Eq.(17), we obtain 

f ^ 

p%^{x) = dhY, Pm{A; b)F^"^'> ix), (19) 

m=l 

where the passing-through term is omitted. Eq.(19) agrees with the result in [12]. 
Similarly, taking the < Ml\ . . .\N > matrix element of Eq.(18), we get 



A-l I A-l-1 



pZi{x)= [dbJ2Yl E A„(6)P/^(i;6)A(^-/-l;6) 



1=0 s=0 t=0 



X yi Jvi ^2 Jy^ 2/3 2/1 2/2 



where ^("^(0;) and F^*\x) are defined in a similar form of Eq.(13). From Eq.(19), 

we get 

= / ^ = / - Xa{b))^ - (1 - A„(^) - Xa{b)n 

(2) N + A ^ N + X and Ml + X for the target fragmentation regions(a; < 0) 

We apply the sum decomposition rule of Eq.(53) in Appendix to Qf^j^{b,a) 
^From Eq.(15), we have 

Spcq1va(&, a) = ^(1 - A„(6))^-iQ„(6, a). (21) 
From Eq.(16), we have 

SpcQ%A{b,a) = -^[1 - (1 - A„(6))^ - A„(6)A(1 - A„(6))^-i]Q„(6, a) 
K[b) 
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+ [1 - (1 - A„(6))^]i(a) + - Xa{b))^-^ + AX^ib) - l]QL{b,a). (22) 

Xa{b) 

Taking the < iV| . . . |A'' > matrix element of Eqs.(21) and (22) except for the 
passing-through terms, we have 

P%^{X) = I db{^\l - (1 - Xa{b))^]F{^) 
J Xaib) 



+ ^[(1 - A„(6))^-i + AX^{b) - l]N{x)}. (23) 
While we take the < M^j . . .\N > matrix element of Eq.(22), we get 

Pj^^(^) = / db [1 - (1 - X,{b))^]A{x). (24) 
From Eq.(24), we get 

^f/ = ^PZiix) = I b[l - (1 - A„)^]. 
Thus, we obtain agi = a^i + = / db[l - (1 - A„(6) - A„(6))^]. 

4 Single-particle distributions in A-A collisions 

4.1 Example(i)-D collisions) 

We consider D-D collisions with A = A ^ 2 in order to show the gross features 
of A-A collisions. All mechanisms in the operator set((a),(b),(c) and (d)) contribute 
to A-A collisions. In the following discussions, we omit the arguments b and a for 
simplicity. We consider the projectile fragmentation regions of D{x > 0). We 
obtain the seven mechanisms as follows; 

(1) Non-annihilation interaction term 



Spr 



Qn Qn 



= 2{l-X,rQi. 



(2) Annihilation terms of order Aa 



Spr 



Qn Qn 



Xa[{l-XafAQn + {l-Xa)Ql], 



(25) 



(26) 



Spr 



a 

Qn Qn 



X4{l-XafQLA + Ql], 



(27) 



Spr 



Ql Qn 
A Qn 



= Xail - Xa)[QLQn + AQn] 



(28) 
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Fig. 4: Diagram for Eq.(31) in D-D collision. The intersection points 
represent M^-M^ collision(x), N-N collision(») and N-N annihilation 
interaction(o). The solid and dashed lines denote nucleon(A'^) and leading 
meson(Mi), respectively. 



Spr 



Qn Ql 

Ql a 



= Xa[QnQL + {l-Xa)QLA]. 



(3) Annihilation terms of order 



Spr 



A Ql 
Ql a 



XI[AQl + QlA], 



Spr 



Mo A 

A Qn 



A2[(l-A„)Moi + iOn], 



(29) 

(30) 
(31) 



where the relations ijl + Xl = ^ and r?m + A^ = 1 are used. The diagram for 
Eq.(31) in D-D collision is shown in Fig.4. 



4.2 Generalization 

We investigate A- A collisions characterized by Q^^(^, a) with A rows and A 
columns. 

(a) Non-annihilation term 

The non-annihilation term in A- A collisions is characterized by Q^ji{b,a). Its 
elements Q'^j with i = 1,2, . . . , A and j — 1,2, . . . , A arc given by 



Then we get 



Q7i =Qn, i = 1, 2, . . . , A and j = 1, 2, . . . , A. 



SprQ^^ = Spr[g^J = A{1 - A„)^ -^(Q„)^. 



(32) 
(33) 



Taking the < ]V| . . . |7V > matrix clement of Eq.(33) except for the passing-through 
term, we get the contribution of the non-annihilation term to the single-particle 
distribution of A + A ^ N + X{x > 0) as follows; 

Ptn^^) = dbAj2 Pm{A^;h)Y,m\A\m,A)F^'\x). (34) 

m=l 1=1 
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This form is similar to Eq.(12) for A + B N + X{x > 0). Prom Eq.(34), the 
non-annihilation cross section is given by 

crli = jj dhY, P^{A^; h) = j db[{l - - (1 - A„ - A„(6))^^]. 

{(3) Annihilation terms of order Xa 

If the annihilation operator A is at the fc-th row and l-th column and its operator 
matrix is given by I), the elements Qlj{k, I) with i = 1, 2, . . . , A and j = 

1, 2, . . . , A of Q^j^{k, I) are given by 

{A, i = k and j = I, 
Ql, i = 1,2, - 1 with j = Z and j = 1,2,...,« - 1 .^g. 
with i = l, ^ ' 

Qni otherwise. 

Then, from Eq.(52) in Appendix, we get 

SprQl^ = K\(k - 1)(1 - A„)-+2(Q„)'-iQL(g„)^-' 

+ (1 - A„)-+'+i(Qi)'-U(Q„)^-' + (A - fc)(l - A„)-+i(Q„)^], (36) 

where uj = A'^ — A — k-l. 

If the annihilation term of order Ao is characterized by Qaaj '^^ S^t the relation 

Qaa = Z]fe=i J2i=i QAAi^^^ 0- 

^Prom the first term and third term in right-hand side of Eq.(36), the contri- 
bution of the annihilation term of order Ao to the single-particle distribution of 
A + A ^ N + X{x > 0) is given by 



xPM - I; b)^^^^{F{a)r{rjL + A^iV(a))(F(a))* + ^(1 - A„)^'-=^^+^ 



X [1 - (1 - A„)^][l + {A- 1)(1 - A„)^ - ^(1 - A„)^-i] ^ Pm{A- b){F{ajr, (37) 

m=l 

where the passing-through term is omitted. 

^From the second term in right-hand side of Eq.(36), the contribution of the 
annihilation term of order Xa to the single-particle distribution of A -|- A — ^ Ml + 
X{x > 0) as follows; 

4lai(^") = EEa„(i - A„)^^-^-'=+i < M^|(Q^)'-M(g„)^-'|7V >, 
fe=i 1=1 



11 



A l-l A- 



(1 - Kr-'^^' [1 - (1 - A.)^] E E E - 1' ^^)^*(^ - 



;=i s=o t=o 



X (L(a))M(a)(F(a))*. 



(38) 



(7) Annihilation terms of order 



a2 



In this case, the annihilation collision of order A^ is characterized by Qaa^)^" ^' ■ 
and Qaa{^^ ^'i 0- The elements of Qff{k, I; k', I') of Qaa{^^ ^i'® given by 



'A, i = k, j = I and i = = 
Ql, i = l,2,...,A;-l withj = /, j = 1,2,...,/-! 

with i = k,i = l,2,...,k'-l with j = Z' (39) 
and J = 1, 2, . . . , Z' — Iwith i = k' , 
, Qn, otherwise, 



where k < k' and / < I'. Then we get 



. o2 



SprQjAik, I; k', I') = Xilik - 1)(1 - A^)""'"'' {QnY-'QLiQnf-^-^QLiQn) 



\A-l' 



+ (1 - Xaf-''+HQLy-^MQny'-'-^QL{Qn)^-'' 



+ {k' -k- 1)(1 - A„) 



n+A-i-i'+i/A 



{Qny-'QLiQnY 



+ (1 - Xaf-'+\QLy'-'A{Qn)^-'' +{A- k'){l - A„)"-''-'+2(g„)^], (40) 



n-i'-i+2/A 



La2 



where = - A - fc - fc'. While, the elements of Qffik, I'; k', I) of QjAi^, I'; k', I) 
are given by 



Qtf{k,l';k',l)=i 



'A, i = k,j = V and i = k' ,j = I, 

-Mo, i = k and j = I, 

Ql, z = 1,2,. . .,fc' - l(i ^ A:) with j = /, i = 1, 2, . . . 

withj=r, j = l,2,...,r(jV0 
with i — k,j = 1,2, ... ,1 — Iwith i = fc', 

, Q„, otherwise, 



(41) 



where k < k' and I < I'. Then we get 



SprQf^(fc, I'; k', I) = A„[(fc - 1)(1 - A„)"-'-''+5(Q„)'-iQ^(g„)''-'-iQi(Q„)^-'' 



+(1 - A„)"-'+3(gi)'-iMo(gL)''-'-ii(Q„)^-'' 



+{k' -k- i)(i - A„)"-'+4(g„)'-igi(Q„) 



(1 - Xaf-''+HQLy-'AiQ„)^-' + iA- fc')(l - A,)"-'-''+3(g„)-4]. (42) 



si-;-/'+3//) \Ai 
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If the annihilation term of order is characterized by Q^^, we get 

A—l A—1 A A 

QfA = Y.Y. E Y.lQ'^AA{k,l;k',l') + QfAk,l';k',l)]. 

k=l 1=1 k' = k+l l'=l + l 

^From Eqs.(40) and (42), wc may calculate the single-particle distributions of 
N and Ml in the projectile fragmentation regions of A by means of the following 
relations: 

'^j|^2(^' ^) ~< -^|SprQ^^|iV > — < passing — through term >, 

Pita2ib,a) =< ML\SpvQf^\N > . 
{S) Annihilation terms of order 

We consider the annihilation collision of order A^. We investigate the single- 
particle distribution of the leading meson Ml- We estimate the following two cases: 

(1) Qaa whose matrix elements Q^f'^ are defined by 

(jaA,! ^ / ^' ^ = 

I Ql, ij^j, 
where i = 1,2, . . . ,A and j = 1,2, . . . ,A. Then we get 

SprQll'' = A^Xi(^^)'^(^^)^"'"'- 

1=0 

(2) Q^jI whose matrix elements Qif'^^ are defined by 

r A, i+j = A + i, 

Q-^'"=<^ Afo, i = l,2,...,A-landj = l,2,...,A-i, 
[ Qn, i = 2,3,...,Aa.ndj = A-i + 2,...,A. 

Then we get 

Taking the < Ml \ . . .\N > matrix element of Eq.(44), we obtain 

1=0 s=0 t=0 

^ f'dm f'dy, ^i(«)(iL)^(yi)iVW(y3). (47) 

Jx yi Jyi Jy2 2/3 Vl 2/2 
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(43) 



(44) 



(45) 



(46) 



While we take the < Ml \ . . .\N > matrix element of Eq.(46), we obtain 

xP*(A-/-l;6) f'^ f'-^M^rn\-)A{'-^)F^Hy.\ (48) 

Jx yi Jyi 2/2 Jy^ VS yi 2/2 

where Pr{l;b) = (^'^ Vm{by-' Xm{by . 

We obtain the similar results for the target fragmentation regions (a; < 0). 

5 Conclusion and Discussion 

We have investigated the single-particle distribution of N, N and Ml in A- 
A and A-A collisions on the basis of EMCM. By introduction of the operator 
matrix with the rule of the sum decomposition, a unified treatment of A-A and 
A-A collisions at high energy is given. The analytic forms of the single-particle 
distribution of the inclusive process A{A) + B ^ C + X are derived. 

Wc have mainly discussed the peripheral reactions. The introduction of the non- 
peripheral reaction Ml + N{N) — >■ N{N) + Ml+ hadrons gives the contributions of 
the higher order of Xa to the inclusive distributions of A-A collisions. Furthermore, 
the inclusive distributions of the secondary particles and the transverse momentum 
distribution are the unsolved problems. Also, the investigation of the relation to 
the Monte Carlo simulation models is one of the interesting problems in future. 



Appendix Definition and properties of operator matrix Oab (ct) 

The operator matrix Oab {oi) in the moment space is defined as the rectangular 
or square array of the operators Oij (a) with A rows and B columns 

dii(a) di2{a) ■■■ OiB(a) 



OAB{a) = 



_dAi{a) dA2{a) ■■■ Oab (a) 



(49) 



by 



The operator matrix OAsict) has the following properties: 

(i) The multiplication of Oij{a) by the scalar quantity (probability) aij is defined 



aiidn{a) 



aAiOAiia) 



aindiBia) 



aABOAB{a) 



A B 

i=lj=l 



On (a) 



OAi{a) 



diB{a) 
Oab (a) 



(50) 



(ii) The rule of the operator addition is defined by 
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Oii(a)+Ai(a) 012(a) ••■ Ois(a) 
021(a) 622(0;) ••• 62b(q;) 

Oai(q;) OA2{a) ■■■ 6^5(0;) 



On (a) 012(a) 
021(a) 622(a) 

6Ai(a) 6^2(0) 



OiB(a) 
02B(a) 

6^B(a) 



Pii(a) 012(a) ••• OiB(a) 
621(a) 622(a) ••• 62B(a) 

6^1 (a) 6A2(a) •■• OABia) 



.(51) 



(iii) The sum decomposition rule of Oab{c() 

(a) The sum of the products of the row elements(SprOAB(a)) 

A B 

(b) The sum of the products of the column elements(SpcOAs(a)) 

B A 

Spc6^B(a)=^[]6i,(a). 

j=i i=i 



(52) 



(53) 



The factorization rule of Eq.(50) is very important. Prom Eq.(50), we may 
calculate the event probability in A-B collisions correctly. 

Notice that the sum decomposition rule must be used after the scalar factors 
taken outside. Namely, for Eq.(50), we have 



Spr 



aiiOii(a) 



aiB6iB(a) 



A B 

i=i j=i 



aAiOAi{a) •■• aABOABia) 
"611(a) •■• 6iB(a)" 



6^1 (a) 



OAB{a) 



For Eq.(51), we have 



(54) 



Spr 



611(a) +Ai(a) 612(a) ••■ 6iB(a) 
621(a) 622(a) ■•■ 62B(a) 

6Ai(a) 6^2(0) ••• 6AB(a) 



Spr 



611(a) 612(a) •■• 6iB(a) 
621(a) 622(a) •■• 62B(a) 

OAi(a) OA2(a) ••• OABia) 
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Spr 



Piiia) 012(a) 
621(a) 622(a) 

OAi(a) 6^2 (a) 



Oii3(a) 
62B(a) 

OAB{a) 



(55) 



Also, Eq.(53) satisfies the similar relations. According to these rules, it is easy 
to show the following relations for Eq.(49) with 6^ (a) — Qo{a) ior i — 1,2, . . . , A 
and j = 1,2, . . . , B where Qo{a) ~ ?7G(a) + A J(a): 



Spr 



60(a) 60(a) 
60(a) 60(a) 



60(a) 
60(a) 



A{r, + X) 



AB~B 



(56) 
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